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Beams usually do not travel through the magnet centers due to errors in storage rings. The beam deviating 
from the quadrupole centers is affected by additional dipole fields due to magnetic field feed-down. The beam- 
based alignment (BBA) is often performed to find a golden orbit, on which the beam circulates around the 
quadrupole center axes. For storage rings with a large number of quadrupoles, the conventional BBA procedure 
is time-consuming, especially in the commissioning phase due to the necessary iterative process. Additionally, 
the conventional BBA method can be affected by strong coupling and nonlinearity of the storage ring optics. 
In this work, a novel method based on a neural network is proposed to find the golden orbit in a much shorter 
time with reasonable accuracy. This golden orbit can be directly used for operation, or can be adopted as the 
starting point for the conventional BBA. The method is demonstrated in the HLS-II storage ring for the first 
time, through simulation and online experiments. The results of the experiments show that the golden orbit 
obtained using this new method is consistent with that from the conventional BBA. The development of this 
new method and corresponding experiments are reported in this paper. 
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I. INTRODUCTION 


Ideally, the beam in a storage ring should circulate on the 
orbit passing through the axes of all magnet centers, which 
is called the golden orbit. The beam orbit may deviate from 
the ideal path due to various errors, such as misalignment, 
magnet imperfection, power regulation errors, etc. When the 
beam traverses the magnets with orbit offsets, it will see unde- 
sired magnetic fields which is called feed-down [1]. The feed- 
down of a quadrupole with an orbit offset causes an additional 
dipole field. To minimize this effect, the beam-based align- 
ment can be adopted to determine the golden orbit for ma- 
chine operation. It has been widely used in the commission- 
ing of storage rings [2, 3]. For storage rings with long circum- 
ference, such as most diffraction-limited storage rings (DL- 
SRs), the number of quadrupoles is large and the conventional 
BBA method becomes much time-consuming [4]. Recently a 
fast BBA method is developed in the ALBA light source using 
AC excitation of the orbit correctors and fast beam position 
data acquisition [5-7]. At HLS-II, with no need to upgrade 
the hardware, a machine learning (ML) based method is de- 
veloped to find the golden orbit for storage rings[8]. 

Neural networks (NNs) have been widely applied in ar- 
tificial intelligence and have gained great success in vari- 
ous fields. Its application has also been introduced to the 
area of particle accelerators [9-12]. At Advanced Light 
Source (ALS), an NN model is used to maintain the verti- 
cal beam size when the gap of insertion devices varies [13]. 
The NN model can also be used to greatly reduce the simu- 
lation time for the optimization of beam dynamics [14]. At 
Shanghai Synchrotron Radiation Facility (SSRF), the image 
processing technique using convolutional neural networks 
(CNNs) is adopted to extract bunch longitudinal phase infor- 
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mation [15]. These applications show great potential of NNs 
in improving accelerator performances. 


In this paper, we present a new BBA method that uses an 
NN model to predict the golden orbit of a storage ring. To 
initiate the experiment, different closed orbits are generated 
by randomly changing the strength of all orbit correctors. 
The beam with various orbit deviations in the quadrupoles 
is subject to various amount of influence from their feed- 
down. This effect can be evaluated by measuring the orbit 
change caused by quadrupole strength variation. Since the 
beam on the ideal golden orbit should not be disturbed by 
changing the quadrupole strength, an NN model can then be 
trained to search for the orbit that is least affected by varying 
quadrupole strength. To train this model, the orbit differences 
due to quadrupole change is used as the input data, and the 
corresponding orbits before quadrupole adjustment is used as 
the output data. The golden orbit is then predicted by setting 
the input value as zero to the NN model. 


The new BBA method is tested in the HLS-II storage ring 
through simulation and online experiments to demonstrate its 
validity. The result shows that the golden orbit obtained from 
the NN model is consistent with that obtained by several it- 
erations with the conventional method. The golden orbit ob- 
tained using this method can be directly used for operation 
or used as a Starting point to speed up the conventional BBA 
that requires several iterations. In general, this new method 
is less time-consuming compared to the conventional BBA 
especially in the case during the initial commissioning [16]. 


In the following sections, the methods of the conventional 
BBA and the new BBA using an NN model are shown in Sec- 
tion II. The simulated result using these two BBA methods in 
the HLS-II storage ring is described in Section III. The on- 
line experiments using both BBA methods are introduced in 
Section IV. Finally the work is summarized in Section V. 
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II. BEAM-BASED ALIGNMENT 


A. Conventional BBA method 


The purpose of the BBA is to find a reference orbit in which 
the beam passes the centers of all quadrupoles in a storage 
ring using beam position monitors (BPMs) and orbit corrector 
magnets (OCMs). A quadrupole generates dipole fields with 
strengths 


Bz = Bopo Koyo, 
By = BopoKoxo, 


() 
(2) 
where Bopo is the magnetic rigidity, Ko is the normalized 
quadrupole strength, x9 and yo are the beam offsets relative 
to the quadrupole center in the horizontal and vertical plane, 


respectively. Therefore, changing the quadrupole strength by 
AK causes a dipole field variation by 


AB, = BopoAK yo, 
AB, = BopoAK 20, 


(3) 
(4) 


resulting in a kick which leads to an orbit change at an obser- 
vation point s by [17] 
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(5) 
VOB 


2sin(7v) 


cos(|¢(s) — ¢(so)| — mv), 


where Lo is the length of the quadrupole, v is the betatron 
tune, 3(s 9) and 3(s) are the beta functions at the locations of 
the quadrupole and observation point respectively, (59) and 
@(s) are the phase advances at the locations of the quadrupole 
and observation point and u stands for the beam positions in 
the horizontal and vertical plane. This equation shows that the 
beam orbit can be affected by the quadrupole strength varia- 
tion and also the beam positions in the quadrupoles. To avoid 
this effect, the reference orbit of the orbit feedback system is 
usually set to the centers of all quadrupoles with u = 0. This 
reference orbit can be determined using the BBA technique. 

The quadrupole center is measured using its nearest BPM. 
Suppose that when the beam goes through the quadrupole 
center, the related reading of this BPM is vo. According to 
Eq. (5), by changing the quadrupole strength AK, the beam 
orbit change is given by 


Au = AKF(v — vo), (6) 
where v is the reading of the target BPM before quadrupole 
strength change and F is the coefficient which can be eas- 
ily obtained from Eq. (5). To measure the quadrupole cen- 
ter, the beam is set to several difference positions at its re- 
lated BPM. For each position, the quadrupole strength is 
then changed with the same AK and the corresponding or- 
bit change is recorded. By applying a linear fit to Eq. (6), the 
quadrupole center vo is obtained. This conventional BBA 
always determines the horizontal and vertical offsets sepa- 
rately [18]. The above analysis implies that the coefficient 
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F is treated as a constant, which means the beam optics re- 
main unchanged during BBA process. In fact, the change 
of quadrupole strength and the closed orbit distortion can af- 
fect the beam optics. At the beginning of commissioning, the 
beam orbit and beam optics are possibly far different from 
the ideal model which induces strong nonlinearity and cou- 
pling. In this case, it needs several iterations for the con- 
ventional BBA method to eliminate the nonlinear effects and 
obtain more accurate quadrupole centers. A neural network 
with multi-layers which deals with nonlinear problems can be 
adopted for the BBA process [19, 20]. 


B. BBA using a neural network 


BBA is based on the principle that the off-axis beam is af- 
fected by the quadrupole strength change. The golden or- 
bit can then be evaluated using the relation between the or- 
bit changes and the initial beam orbits before varying the 
quadrupole strength. This relation can be explored by train- 
ing a neural network using the orbit changes as the input data 
and the initial orbits as the output data. By setting the orbit 
change to zero, the corresponding initial beam orbit is just 
the predicted golden orbit. The main idea of this new BBA 
method is shown in Fig. |. To obtain data for training the NN 
model, the simulation or online experiment is carried out as 
follows: 


e Randomly exciting all corrector magnets to form a ini- 
tial closed orbit; 


e Recording all BPM readings; 


e Changing all quadrupoles with a same amount to form 
a new closed orbit. 


e Recording the changes in all BPM readings; 


e Resuming the quadrupole and corrector strengths to the 
original values; 


e Repeating the above procedures. 


A typical dense neural network has one input layer, sev- 
eral hidden layers (also called middle layers), and one output 
layer, as shown in Fig. 2 [21]. The nodes are called neurons 
where the data are transferred. The nodes between adjacent 
layers are connected to each other by an arrow, which shows 
the flow of data. Each arrow represents a linear transforma- 
tion combined with an activation function used to introduce 
nonlinearity if necessary [11]. A loss function is used to de- 
scribe the performance of the neural network. The NN also 
needs an optimizer function to optimize the parameters used 
for data transmission. The optimization is carried out by min- 
imizing the loss value. 


Il. SIMULATION STUDY FOR THE HLS-II STORAGE 


RING 


Simulation is carried out to evaluate the validity of the new 
BBA method based on an NN model before online experi- 
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Fig. 1. Schematic of the neural network based BBA method. Differ- 
ent orbits are generated by randomly adjusting the orbit correctors. 
On each orbit, all quadrupoles are changed with the same AK at 
the same time to generate orbit changes. The orbit changes are used 
as the input data of the neural network and the corresponding initial 
orbits are used as the output data for the training. 
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Fig. 2. Diagram of a typical dense neural network, which consists 
of one input layer, several hidden layers and one output layer. Here 
hyperbolic tangent (tanh) is adopted as the activation function. 


ments. The accelerator toolbox (AT) is used for the simula- 
tion in this work [22]. TensorFlow, which is adopted in this 
work, provides a flexible platform that makes it easy for users 
to build and train an NN model [23, 24]. 


The HLS-II storage ring has two super periods with a cir- 
cumference of 66.1m. The layout for one super period is 
shown in Fig. 3. The orbit system of the storage ring consists 
of 32 BPMs and 32 correctors combined to the sextupoles. 
32 quadrupoles are installed which need to measure their real 
centers through the BBA procedure [25]. 


Random rotation and shift errors are applied to simulate the 
misalignment of the elements and girders. The errors are gen- 
erated in a normal distribution with truncation at three stan- 
dard deviations. According to the design report, the error set- 
tings of all magnets, girders and BPMs are listed in Table 1. 
A set of misalignment errors of the whole ring is shown in Ta- 
ble 4. The magnet strength errors are also applied. The BPM 
random measurement error is set as 0.5 um [26]. 


166 


188 


B functions (m) 
Horizontal 7 function (m) 


Fig. 3. One super period of the HLS-II lattice. There are 32 
quadrupoles and 32 BPMs in the storage ring. The 32 combined- 
function sextupoles are used as the horizontal and vertical correc- 
tors. 


Table 1. Misalignment error settings used for the HLS-II storage 
ring. 


Type Shift error (um) Rotation error (urad) 
Axis X Y S X Y S 
Girder 50 50 200 500 500 500 
Dipole 200 200 150 500 500 500 
Quadrupole 200 200 150 500 500 500 
Sextupole 200 200 150 500 500 500 
BPM 200 200 150 500 500 500 


A. Conventional BBA method 


In the conventional BBA measurement for one quadrupole, 
the beam is moved to three different positions with the help 
of the corrector magnets [27]. At each position, the change 
in the beam orbit from all BPM readings is recorded after 
varying the strength of the target quadrupole with a certain 
AK. The orbit changes can be fitted linearly as a function 
of the beam position in the target quadrupole. An immobile 
point can be found by setting the position at which the BPM 
changes vanish. The quadrupole center is then obtained by 
adding up all immobile points from each BPM. A whole BBA 
routine repeats this process for all quadrupoles in both hori- 
zontal and vertical planes in the storage ring. To increase the 
BBA accuracy, one can repeat the measurement after moving 
the beam to the orbit obtained from the previous BBA exper- 
iment. This scheme is usually needed at the machine com- 
missioning stage. Fig. 5 shows the simulated measurement of 
the horizontal and vertical centers of one quadrupole magnet 
in the HLS-II storage ring. At least three conventional BBA 
iterations are needed to decrease the standard deviation of the 
fitted Gaussion function of the quadrupole center to several 
microns, which is the same order of BPM measurement reso- 
lution [28]. 
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Fig. 4. Misalignment errors applied to the HLS-II storage ring. (a) 
Shift errors. (b) Rotation errors. 


B. BBA using an NN model 


In the simulation, the correctors are set randomly with a 
certain range kicks to move the beam orbit. Here, the kick 
angle variations are generated by a normal distribution with 
the standard deviation of 0.05 mrad and a truncation at three 
standard deviations is applied. For each random orbit, all 
quadrupoles are simultaneously changed with a same amount 
of AK (—0.02m~?). The corresponding initial beam orbit 
and orbit changes are recorded from all BPMs. 

The whole simulation generates 10000 samples. In each 
sample, there are 64 initial orbits and 64 orbit change data, 
with 32 in the horizontal plane and 32 in the vertical plane. 
These samples are adopted for training the neural network. 
Fig. 6 shows the random initial beam orbits within a range 
of (—5,5) mm. Fig. 7 shows the orbit change after the 
quadrupole adjustment. The range of orbit change is within 
(—1.5,1.3) mm and (—0.8,0.8) mm in the horizontal and 
vertical plane, respectively. 


To obtain the golden orbit, an NN model is trained using 
these data. The 64 sets of orbit change data are set as the 
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Fig. 5. Simulated BBA measurement in the HLS-II storage ring. 
The horizontal and vertical measurement is applied respectively. (a) 
Horizontal quadrupole center measurement. (b) Vertical qaudrupole 
center measurement. The plots show the orbit change observed from 
all BPMs by varying the target quadrupole strength with a certain 
AK when the beam is at three different positions. For each BPM, 
its changes can be fitted to find a center for the quadrupole. All of 
the found centers are then fitted using the Gaussian function. The 
red line shows the fitted centers using all BPMs, which represents 
the BBA center of this quadrupole. 


input of the model, and the 64 sets of the corresponding ini- 
tial orbit data are set as the output of the model. 80% of the 
data are used to train the model and the rest are used to test 
the performance of the model. There are 128, 256 and 128 
neurons in the three hidden layer, respectively. The tanh is 
taken as the activation function to provide nonlinearity. The 
NN model is trained using the Adam optimizer [29]. The loss 
function is the mean squared error (MSE) between the mea- 
sured data and model-predicted result, which is 


(7) 


Fig. 8 shows the comparison of the golden orbits obtained 
from the conventional BBA and the BBA using a neural net- 
work. The good consistency of these two BBA methods 


loss = mean((Fmodel — Yreal)”)- 
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Fig. 6. Distribution of the initial orbits generated by randomly 
adjusting the orbit correctors within a certain range. (a) Initial orbits 
in the horizontal plane. (b) Initial orbits in the vertical plane. 


shows the validity and effectiveness of this new BBA tech- 
nique. Next the online experiment is carried out in the HLS-II 
storage ring. 


IV. ONLINE EXPERIMENT IN THE HLS-I STORAGE 
RING 


The conventional BBA has been applied to the HLS-II 
storage ring [30]. Fig. 9 shows the BBA result for one 
quadrupole. The errors of fitting for most quadrupoles are 
within 20 um. 


A. Training data acquisition 


Similar to the simulation, the training data could be ob- 
tained from the real storage ring. Before experiment, the 
magnet strengths are set according to the result from the early 
commissioning of the storage ring. In this case, the beam 
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Fig. 7. Distribution of orbit change caused by varying the 
quadrupole strength for each random orbit in the simulation. (a) 
Horizontal orbit changes. (b) Vertical orbit changes. 
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BBA and the NN-based BBA. The simulation shows good consis- 
tency of these two methods. 
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storage ring. The upper and lower plots show the horizontal and 
vertical orbit changes observed from all BPMs by varying the target 
quadrupole strength with a certain AK when the beam is at three 
different positions. For each BPM, the change in its reading can be 
fitted to find a center for the quadrupole. The red line shows the 
averaged fitted centers using all BPMs, which represents the BBA 
center of this quadrupole. 


is not on the golden orbit which connects the centers of the 
quadrupoles. The online experiment for obtaining training 
data is reported in this subsection. 


During the online experiment, the orbit feedback system 
is turned off and the correctors are randomly set to gener- 
ate different orbits. As a compromise of beam stability and 
data diversity, the adjustment range of all correctors is set to 
+0.8 A relative to the starting point. This adjustment range 
of correctors ensures no beam loss during the experiment by 
controlling the orbit change within a distinguishable range, as 
shown in Fig. 10. The horizontal tune in the HLS-II storage 
ring is about 4.44 while the vertical tune is 2.36 which is fur- 
ther from the half-integer. When the quadrupole strength is 
increased simultaneously, the horizontal tune increases and 
easily reaches the half-integer resonance and causes beam 
loss. Therefore, all quadrupoles strength is adjusted in the 
decreasing direction, by an amount of —0.02 m~? (normal- 
ized focusing strength). After the orbit change is recorded, all 
quadrupole strength is restored to the original values. For the 
HLS-II storage ring, the time constant for the orbit corrector 
power supplies is about 15 ms [30]. The time constant for the 
quadrupole power supplies is about 30 ms. This means that 
one complete loop of this measurement could be done within 
1 second. To ensure the accuracy of the data acquisition, the 
measurement time for one loop is set to 2 seconds. 


The online experiment is carried out during the machine 
study time [31]. The whole measurement generates 21000 
samples. These samples are adopted for training the neural 
network. Fig. 10 shows the randomly generated initial beam 
orbits before varying the quadrupole strength. The distribu- 
tion shows that the orbits are generated within the range of 
about (-10,10) mm, and the densest distribution is around 0. 
The orbit change after the quadrupole adjustment is also an- 


270 alyzed, and the distribution of the orbit differences is plotted 
27 in Fig. 11. The range of orbit change is within (-3,2) mm in 
272 the horizontal plane and within (-1.5,1.5) mm in the vertical 
273 plane. 
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Fig. 10. Distribution of the orbits generated by randomly adjusting 
the orbit correctors within a certain range. (a) Horizontal BPM read- 
ings. (b) Vertical BPM readings. 


B. NN model training using online data 


In this subsection, the relation between the initial orbit and 
orbit change after quadrupole adjustment is explored using a 
dense neural network. Similar to the simulation, the 64 sets 
of the orbit change data are set as the input to the model, and 
the 64 sets of the corresponding initial orbit data are set as the 
output of the model. To figure out the data size requirement, 
two models are trained with different numbers of samples. In 
model I, all 21000 samples are adopted, 5/6 of the samples 
are used for training, and 1/6 are the validation data set. As 
a comparison, model II is trained with only 3000 samples in 
the training set and 600 samples in the validation set, namely 
3600 samples are adopted in total. The Adam optimizer and 
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Fig. 11. Distribution of orbit change after varying the quadrupoles 
for each random orbit. (a) Horizontal orbit change. (b) Vertical orbit 
change. 


the MSE loss function are used for training models. 

The trained NN models are evaluated by calculating the 
mean absolute error (MAE) between the measured real values 
and model-predicted values of the validation samples for each 
BPM: 


(8) 


The absolute error for both models is shown in Fig. 12, which 
shows that the errors of model I are smaller than those of 
model II. In the horizontal plane, the overall averaged abso- 
lute error is around 99 um for model I and 125 um for model 
II. In the vertical plane, the overall averaged absolute error is 
around 62 um for model I and 71 um for model II. The result 
shows that increasing samples for the NN model training can 
improve the model accuracy. 


MAE = mean (|T measured = F predicted] )- 


C. Golden orbit from the NN model 


In the NN training, the orbit changes caused by varying 
the quadrupoles are used as the input data. The correspond- 
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Fig. 12. The mean absolute error (MAE) between the measured 
beam orbits and predicted values of the validation samples for all 
BPMs. (a) Horizontal plane. (b) Vertical plane. 


ing initial orbits are used as the output data. The beam on 
the golden orbit should have the least orbit distortion (ideally 
zero) due to the change in quadrupole strength. Therefore, we 
can set the input as zero to the NN model and the correspond- 
ing output is just the golden orbit. 


To estimate its accuracy, this golden orbit is compared with 
that obtained using the conventional BBA method and the re- 
sult is shown in Fig. 13. The sub-figures in Fig. 13 show 
the difference between the novel and conventional BBA. The 
result shows that this golden orbit is consistent with that ob- 
tained from the conventional BBA. In the horizontal plane, 
the averaged difference between the conventional BBA and 
the model prediction is around 46 um for model I and 53 um 
for model II. In the vertical plane, the averaged difference 
between the conventional BBA and the model prediction is 
around 42 um for model I and 39 um for model II. 


Although the training error of model I is obviously smaller 
than that of model II, the difference in the predicted golden 
orbits from these two models does not have a large devia- 
tion [32]. In the HLS-II storage ring, the typical experimental 
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Fig. 13. Comparison of the golden orbit obtained using the NN 
model and the conventional BBA method. The difference between 
the two golden orbits is shown in the subfigure. (a) Horizontal 
golden orbit. (b) Vertical golden orbit. 


period for the conventional BBA process is around 5 hours. In 
the machine commissioning phase, this BBA process needs 
to be repeated several times to obtain a precise result. Model 
II uses only 3600 samples, which leads to a shorter online 
measurement time (about 2 hours). As discussed above, the 
online measurement time for this new method is irrelevant to 
the total quadrupole number. This differs from the conven- 
tional BBA where the larger the storage ring, the more time it 
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takes. On the other hand, the NN trained golden orbit can be 
set as the starting point for the conventional BBA. This helps 
to reduce the iterative process of the BBA starting from the 
initial commissioning orbit and hence the experimental time. 


V. SUMMARY 


A novel method is developed to search for the golden or- 
bit of a storage ring. This method trains a neural network 
model using the simulated or online data of different closed 
orbits and the corresponding orbit change caused by varying 
all quadrupole strength simultaneously. The online experi- 
ment can be done with less time especially for large storage 
rings. This golden orbit is compared with the one obtained 
using conventional BBA and the result shows good consis- 
tency. 

The NN-based BBA is a good choice for the commission- 
ing stage of a storage ring where the beam optics is far differ- 
ent from the ideal model and the closed orbit is deviated from 
the magnet centers. In this case, the linear process of the con- 
ventional BBA is not accurate anymore. Besides, the conven- 
tional BBA treats the horizontal and vertical orbits separately. 
However, the coupling of a real machine is non-negligible es- 
pecially when the coupling is not sufficiently corrected. The 
NN-based method deals with the transverse planes simulta- 
neously which naturally solves the coupling issue. Addition- 
ally, the new BBA method can be better applied to storage 
rings with strong nonlinear effects, which is often the case 
for DLSRs. With strong nonlinearity, the conventional BBA 
method might work within a limited region since the linear- 
ity of the orbit response is assumed. Since the NNs can be 
used to solve nonlinear problems as is well known, this NN- 
based BBA method is expected to work more effectively for 
the DLSRs. In another perspective, this new technique can 
better deal with the cases where the quadrupoles are pow- 
ered in series, since there is no need to vary the strength of 
all quadrupoles individually. It is anticipated that some small 
light sources or boosters will be benefited from this new BBA 
method. 
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